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Abstract: 

A computer program for evaluating colour factors of QCD Feynman diagrams is presented, 
and illustrative examples on how to use the program to calculate non trivial colour factors are 
given. The program and the discussion in this paper is based on a diagrammatic approach 
to colour factors. 
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Program Summary 



Title of program: colour 

Program obtainable from: |http: / / www.thep.lu.se/tf2 /hep 

or ftp:thep.lu.Se/pub/LundPrograms/Misc/colour.tar.Z 

Computer for which program is designed: Computers with an ANSI C compiler. 

Program language used: C 

Peripherals used: terminal and mass storage for input, terminal for output 
Number of lines in combined program: ^1000 
Keywords: Feynman diagrams, QCD colour factors 

Nature of the physical problem: In many QCD Feynman diagrams there is a need to calculate 
the colour factor. This program calculates colour factors fast and accurately. 

Method of solution: QCD Feynman diagrams, with no four-gluon vertices, factorize into a 
colour (group) factor and a kinematical factor. The colour part of any closed QCD Feynman 
diagram can be transformed into a sum of diagrams including only closed quark (colour) 
lines. In each such term, the number of quark (colour) loops is counted, giving factors of 3 
(or N c ). The sum of these terms is then the desired colour factor. 

To perform the translation, we need to have rules to interpret gluons and vertices into quark 
lines. There is in principle only need for three equalities to be able to calculate any closed 
QCD diagram (not containing four gluon vertices). To achieve good computing performance 
further equalities have been introduced. The most essential ones are listed in Appendix [A|. 

Restriction of complexity of the problem: Four gluon vertices cannot be included directly in 
diagrams to be calculated. The number of gluons in diagrams to be calculated are limited 
to 200, but can easily be changed prior to compiling. 

Typical running time: For typical diagrams, fractions of a second. Complicated diagrams - 
a few seconds to hours. 
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1 Introduction 



There are several programs available to perform momentum Feynman diagram calculations, 
either symbolically or numerically JTJ. Some of these perform QCD calculations, but usually 
the QCD calculating programs do not calculate the colour factors for a given diagram. The 
users have to calculate these factors themselves. This can be done in different ways, the 
most appealing method is presented in ||, Q where only diagrammatic manipulations are 
needed to calculate colour factors of diagrams. 

In this paper we present a computer program for calculating colour factors of QCD diagrams 
using the rules and equalities derived in [Q. An introduction is given on how to use the 
program and which results can be obtained. Using a set of projection operators non trivial 
results can be achieved, such as how the interaction probability between two gluons depend 
on their common colour state. Only closed Feynman diagrams (vacuum QCD diagrams, ie. 
diagrams with no external legs) can be calculated, but this is not a real problem since the 
restriction can be circumvented as will be shown. We also present a short introduction to 
the group theory part of the calculations. 



2 Theory in brief 

The program presented here is based on the results in and ||. These results has been 
further developed and explained to us by Yuri Dokshitzer and Bo Soderberg [ffl. In this 
section and the appendices we give a brief presentation of the main results. 

The existence of colour factors in QCD originates from the invariance of the Lagrangian un- 
der colour transformations. For a certain interaction, with incoming and outgoing particles, 
intermediate propagators and different kinds of vertices, one has to sum over all possible 
colours, if they are not explicitly specified. Since the kinematical factor in the interaction 
probability is equal for all these diagrams (excluding diagrams containing four-gluon ver- 
tices), we therefore have a multiplicative factor to the probability, which we get by counting 
the colours. 

Two preliminary remarks. When we discuss colours and anti-colours we often use quarks 
and anti-quarks, respectively, to describe the colour flow. In most expressions we keep N c as 
a free variable in order to be general. Whenever the number of colours is specified, N c = 3 
has been used. 

Every 3 m (8> 3 n space can be divided into a set of irreducible subspaces. In the simplest 
example of a quark and an anti-quark system, the irreducible representation is the singlet 
and the octet state, ie. 3 <E> 3 = 1 © 8, where the octet state corresponds to the gluons. The 
equation 
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states this fact in equation form. (The factor of 2 in the second term is due to the Gell- 
Mann normalization of gluons, see discussion in Appendix |B[ We have adopted Gell-Mann 
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normalization in all calculations in this paper.) In the LHS of Eq (|T|) a colour and an anti- 
colour state, 3 <S> 3, are divided into irreducible singlet and octet states on the RHS. The 
degeneracy (multiplicity) of a state is given by taking the trace of the representation of the 
state. Calculating the trace in this diagrammatic representation corresponds to connecting 
the incoming legs and the outgoing colour flow. For the LHS we find 

Tr [LHS] = Tr [ ] = = 

as expected, and for the RHS 



Tr [RHS] = 




= 1 + (N*- 1) 



Another observation is that the two diagrams, Pi and P2, in the RHS of Eq (|l|) are disjoint 
projection operators since the terms satisfy Pf = Pi and P{Pj = 0,PjPi = 0. These can be 
used to split any quark-anti-quark diagram into an irreducible representation 

The projection operators can be used to ask specific questions about the colour states par- 
ticles are in, eg. what is the probability to find two free gluons in an anti-symmetric octet 
state. This is easily solved by calculating colour states. The two gluons are in a 8 <8> 8 state 
with degeneracy 64, and the irreducible representation of the 8 <S> 8 space contains an anti- 
symmetric octet (with degeneracy 8), ie. the answer is 8/64 since the probability is equal for 
all states. In the next section we will give a few other examples on these issues and describe 
how to perform the calculations using the program. 

Given a problem concerning colour states, the problem of finding solutions boils down to 
finding the projection operators for reduced representations. Once this is done all interesting 
numbers can be calculated using the program. In Appendix [FJ we have listed the (irreducible 
representation) projection operators for a few states. 

The triple gluon coupling is anti-symmetric in colour. We have generalized the program to 
include also a symmetric coupling between three colour octets ("gluons"). 



3 Using the program 

A straight forward use of the program is to calculate QCD vacuum diagrams. This might 
not be too exciting. Using projection operators more interesting results can be obtained. A 
few illustrative examples will be given below. 

The program interacts with the user through input files. The diagram to be calculated is 
defined in one input file which is read by the program, the diagram is calculated and the 
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result is written on the screen. No interaction with the user is needed after the program has 
been started. 

The command line syntax to start the program is: 

colour [-df file | -h | -help] 

where 

-df Changes the input file name. The default name is colour. in. 

-h or -help Prints a small description of the program and a list of valid options. 



There are few rules for defining a qcd diagram in the input file. 

• Every gluon line is assigned a unique number. 

• Triple anti-symmetric gluon vertices are represented by a letter / followed by the three 
incoming gluon numbers. 

• Symmetric triple vertices^ are represented in the same way as the anti-symmetric 
vertices but with the letter d. 

• Quark loops are represented by a letter t, followed by all gluons connected to the quark 
loop. 

• Four-gluon vertices cannot be directly included in the program since the colour and 
the kinematical parts of the cross section do not factorize. This can, however, be taken 
care of by separating the four-gluon vertex in factorizable terms, as discussed further 
in Appendix |A|. 

• Ghosts can also be included since these couple to gluons with the same colour factor 
as gluons themselves. Thus ghost vertices can be represented with a letter /. 

• For the quark loops and the anti-symmetric gluon vertex the order of the gluons are 
important since f a b c = —fbac- So when translating a diagram into a file all loops and 
/ vertices must be written down in the same order (clockwise or anti-clockwise). 

• Composite objects, d, f, and t, must be separated with a newline character, while gluon 
numbers are separated with a space. All other separation characters are invalid and 
will stop the execution of the program. 



For example, the diagram 



would transform to a file like 

d 1 2 3 
f 2 4 5 
f 3 5 6 
t 1 4 6 

2 As mentioned above these are not physical QCD gluon vertices, but are included for generality. The 
vertex is also used in projection operators. 
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Figure V. Diagrams used to calculate Cp. 



The numbers assigned to every gluon are arbitrary and occur twice in the file since every 
gluon is connected to two vertices in the diagram. If you run this file, the answer will be: 
^3 (iV c 4 - 5iV c 2 + 4) and is equal to 5 using N c = 3. 



3.1 Example 1: Colour factor for gluon emission from quarks 

Let us start with the problem of calculating the colour factor, often denoted Cf, for emission 
of a gluon from a quark line. This can be easily done by hand, but it is a nice example to 
demonstrate how the program can be used. 

We will study a transition matrix T, defined in Fig [l]. A colour singlet comes in, a quark- 
anti-quark pair is created and a gluon is emitted from one of the quark legs. We have chosen 
to use a singlet state since the program can only handle closed diagrams^. By definition, Cf 
is a multiplicative factor to the transition probability, due to the emission of a gluon from the 
quark leg. For the transition probability, we have to square the matrix by multiplying it with 
its conjugate (reversing the colour flow and using the mirror image). The result then must 
be normalized by dividing by the squared amplitude of just creating the quark-anti-quark 
pair. 

Consequently, to obtain Cp we have to run the program twice, one time with the diagram 
TT^ and once with the diagram TqTq. 

In order to use the program we have to translate the diagrams into an input file. This can 
be done as 



Diagram 


Input file 


Result 




t 1 1 






t 





Thus, we get that 



3 There are obviously other diagrams which would enable us to calculate Cf, but these involve obscuring 
technicalities. 
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Pi 





Tt[P,T] = C \ Tr[P] = 



Figure 2: Diagrams used to calculate the singlet eigenvalue, \±, for the state T. 

3.2 Example 2: Partitioning into irreducible representations 

Any diagram in the 8 <8> 8 space can be written as a sum of the projection operators de- 
scribed in Appendix [B] (and crossing operators which are built from operators of the same 
multiplicity). Consider eg the gluon-gluon scattering matrix T in Fig |2|. This matrix can be 
written as a sum of projection operators solely, that is without any crossing operators, 

T = £A^. 

The eigenvalues can easily be found. To be concrete we will determine Ai corresponding to 
the singlet state, and the other eigenvalues are derived in the same way. To calculate this, 
we apply the operator to the matrix T, thus projecting out the singlet part of the diagram. 
Then we take the trace of PT and divide the result with the trace of just the operator, 

TP - A P A - Tr(PlT) 



The translation into a file looks like 



Diagram 


Input file 


Result 


Tr (P X T) 


f 1 2 3 
f 1 3 2 


N c [Nl - 1] 


Tr(Px) 


See Sec |3l] 





Running the program we find that Ai = N c .f\ 



3.3 Example 3: Probability of an octet state 

The question here is: If we have a transition matrix of some kind containing two gluons, what 
is the probability for them to be in any of the states described in Fig |6| in Appendix [FJ. Take 
for example the transition matrix T in Fig [| what is the probability for the two incoming 
gluons to be in the anti-symmetric octet state P 8a ? 

4 Note that the eigenvalues depend on the normalization of the gluon, see Appendix H. 
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T = 




TT 



Tr[rT + ] = 




R = 



1 



8a N„ 




Tr[p R rr + ] = 




Figure 3: Diagrams used in Example 3.3. 



To solve this problem, the first step is to square the matrix, thus getting the symmetric 
matrix TT^ which correspond to the probability. Then the antisymmetric octet part of 
the amplitude is projected out, by applying the operator. The probability is given by the 
multiplicity of the projected part, divided by the total multiplicity of the amplitude, 



Prob(P 8a \T) 



Tr(P 8a TT*) 
Tr(TTt) 



For this we have to run the program twice, once with the diagram Tr(TT') and once with 
the diagram Tr(P 8a TT^). 

The translation into a file looks like 



Diagram 


Input file 


Result 


l 


f 1 4 2 






f 4 7 5 






f 7 9 8 






f 6 8 9 






f 3 5 6 




9 


f 1 2 3 


Nl [AT? - 1] 


iV c Tr(P 8a TTt) 


f 1 4 2 
f 4 7 5 
f 7 9 8 
f 6 8 12 
f 3 5 6 
f 11 2 3 
f 1 10 9 
f 10 11 12 


|[iV c 6 + lliV c 4 -12iV c 2 ] 



Putting it all together gives a probability of 7/24 (with iV c 
gluons in an octet state. 



3) to find the two incoming 
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3.4 Limitations 



There are some limits in the program. First, the physics limitations. Four-gluon vertices are 
not included and therefore diagrams containing these vertices cannot be calculated directly. 
The four-gluon vertex must be replaced with three other diagrams, which have to be cal- 
culated separately (see Appendix There is also one specific diagram that the program 
cannot evaluate, a gluon ring, but the colour factor for this diagram is simply the number 
of different gluons, N% — 1. This is purely technical and the diagram could in principle be 
implemented into the program as a special case, but the gain is not worth the effort. 
Secondly, local implementation limits. There is a limit on the number of gluons which can 
be included in the program. This is set to 200 but can easily be changed by the user before 
compiling. There is also a limit to which power in N c calculations can be done, the default 
range is iV" 100 - N® 9 , this can also easily be changed prior to compiling the program. A 
warning message will appear if over- or underflow occurs. 
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A Colour Rules 



The set of rules in Fig |4], taken from p| , are implemented in the program with the exception 
of the gluon loop Clb. The equation numbers are the ones to be found in 0, and as indicated 
we do not list all of them. We note that in principle it is enough to use equations Cll and 
C19 to translate any diagram without four-gluon vertices to diagrams containing only quark 
loops. The other equalities are implemented in order to increase the computing performance. 



= N c 





N r 2 -1 








J_ 
= C F 



+ 2 



1 

2N C 



2 



2 



2 



Cla 
Clb 
C2 
C7 

Cll 
C14 
C15 

C18 
C19a 



C19b 



Figure 4: Rules for diagram manipulations. 



The four-gluon vertex has not been included in the program since QCD diagrams containing 
four-gluon vertices do not factorize into a colour (group theoretical) part and a kinematical 
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part. The four-gluon vertex, 
has the form 

f ab nfn Cd {g ai g ps - g aS g M ] + fa nd hcn[g ai g ps - g ap g 1& ] + /^W^V 7 - <? a V1- 

Thus a diagram with such a vertex can be split into three pieces, each of which has a factor- 
ized form (this separation is not unique). Each piece has a colour factor corresponding to a 
diagram where the four-gluon vertex is replaced by two three-gluon vertices, as symbolically 
shown in Fig [5[ 

a y^ 8 ^> — £ [g a v 8 - g a v^] + 

j£ [ g «V 5 - g a V 5 ] + 
[g°V Y - g af V 8 ] 




Figure 5: Symbolic representation of the colour factors in a four-gluon vertex. 



B Projection Operators 

In Figs ^-^| we present the set of projection operators needed to obtain irreducible represen- 
tations of the 8 <S> 8, 3 <S> 3, 3 <8> 3, and 3 <8> 3 <8> 3 spaces, respectively. For the derivation of 
these operators we refer to Cvitanovic 0. 

There are different normalizations adopted for gluons and the coupling constant a s . A 
difference in the overall normalization of the projection operators is not a problem as long 
as it is used consistently. The overall normalization is factorized out in most cases, as in 
the Examples |3.2| and |3.3| . In Example |3.1| we have adopted the Gell-Mann (conventional) 
normalization of the gluons. 

The difference in normalization, lies in how we interpret the curly lines, representing the 
gluons, and also in the interpretation of the triple vertices. With the Gell-Mann normaliza- 
tion, the gluon is interpreted as the octet projection in the 3® 3 space together with a factor 
1/2 (compare Eq Cll of Appendix 0). This factor originates from the conventional QCD 
normalization of the generators of the SU(3) group, ie. 

e = (2) 
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A more convenient way to normalize the gluon is to let the factor of | in Eq |2| go into a s , 
ie. to let the gluon represent just the octet projection, without any factors. In the same 
spirit, the triple gluon vertices (symmetric and anti-symmetric) can be defined without any 
factors. 

The program output comes with both normalizations. The answer which is expressed as 
a polynomial in N c , is normalized with the second convention. If you use the polynomial 
ignoring the overall factor (and the numerical answer), you should replace the factors within 
parenthesis, in the operators in Figures || and 0, by a factor 1. If you prefer to use the 
Gell-Mann normalization, you must make use of the factors within parenthesis. 

The projection operator, Pq, has the multiplicity when the number of colours is 3. This 
subspace exists eg. in SU(2) where it has multiplicity 3 and corresponds to a real state. 




Figure 6: Projection operators needed to obtain the irreducible representations of the 
8©8 = 0ffil©8 a ©8 s ffil0©l0©27 space. 



11 



p = — 

1 N c 



P = 



^ (2) 



Figure 7: Projection operators needed to obtain the irreducible representations of the 
3 <S> 3 = 1 © 8 space. 



P = 



P = 



[= x] 

[=-X] 



Figure 8: Projection operators needed to obtain the irreducible representations of the 
3 <S> 3 = 3 © 6 space. 



Pi = 
[ P = 

2 £ = 



6 



6 



io 6 



Figure 9: Projection operators needed to obtain the irreducible representations of the 
3©3©3 = 1©8©8©10 space. 
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